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ABSTRACT 

Numerical simulations of the dynamical friction suffered by a Galactic center star cluster harboring an 
intermediate-mass black hole (IMBH) have been performed. Gerhard has suggested that dynamical friction, 
which causes a cluster to lose orbital energy and spiral in toward the Galactic center, may explain the pres- 
ence of a cluster of very young stars in the central parsec, where star formation might be prohibitively difficult 
because of strong tidal forces. However, numerical simulations by Kim & Morris showed that this is only 
possible if the cluster initially has an extremely dense core. Hansen & Milosavljevic recently suggested that 
the presence of an IMBH in the cluster core might stabilize the core against tidal disruption during the inspiral 
through dynamical friction, and thus might easily deliver young stars down to the central parsec. We find that 
the presence of an IMBH does lower the minimum initial core density required to transport young stars down 
to the central parsec, but this is possible only when the mass of the IMBH is at least ~ 10 % of the total cluster 
mass. This fraction is significantly higher than that estimated by Portegies Zwart & McMillan with numerical 
simulations of IMBH formation by successive merging of stars in the cluster core, so it does not appear that a 
realistic IMBH can help transport young stars into the central parsec. 

Subject headings: stellar dynamics — Galaxy: center — Galaxy: kinematics and dynamics — galaxies: star 
clusters — methods: N-body simulations 



1. INTRODUCTION 

The central parsec of the Galaxy contains a cluster of very 
young stars, including 16 very luminous He I emission line 
stars (Krabbe et al. 1995; Paumard et al. 2001) as well as 
many O and B stai's (Eckait, Ott, & Genzel 1999). Krabbe 
et al. (1995) find that the properties of these stars can be ac- 
counted for by a burst of star formation between 3 and 7 Myr 
ago. The He I stars appear to be evolved massive (> 40 M©) 
stars with stellar ages of ~ 5 Myr (Paumard et al. 2001). 
Krabbe et al. (1995) estimate that the number of OB stars 
with L > 3 X IO^Lq and Wolf-Ray et stars (WNL, WCL, & 
He I stars) in the central parsec is ^ 50. We define these 
objects as Young Massive Stars (YMSs), and estimate their 
progenitor masses to be > 40Mq. Despite their very young 
ages, in situ formation of these stars may be inhibited by the 
strong tidal forces in the central parsec. It is not clear that 
the maximum density of gas currently in the central parsec 
can be as high as the Roche density required for a cloud to 
remain bound, ~ 4 x 10^ Hcm"-'(1 pc/R)^ for galactocen- 
tric radius R. Jackson et al. (1993) report a density of a few 
times 10^ cm"^ at a distance of ^ 1 pc, while Christopher et 
al. (2004) argue that densities approaching the limiting Roche 
density can be found at a distance of ~ 2 pc; (see the discus- 
sion in §|3lbelow) 

One possibility is that the gravitational collapse leading to 
the formation of the present cluster of young stars in the cen- 
tral parsec was triggered by infall of a particularly dense gas 
cloud, which experienced compression by shocks involving 
cloud-cloud collisions, self-intersecting gas streams, or vio- 
lent explosions near or at the central black hole (Morris 1993; 
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Sanders 1998; Morris, Ghez & Beckhn 1999). Alternatively, 
the star cluster could have formed outside the central parsec, 
where tidal forces are relatively weaker and star formation 
is consequently less problematic, and later migrated into the 
Galactic center (GC). Gerhard (2001) proposes that dynam- 
ical friction can bring a massive young star cluster, initially 
embedded in its parent molecular cloud, into the central par- 
sec during the lifetime of its most massive stars, depending 
on the initial location and mass of the cluster. The drag force 
represented by dynamical friction, acting in the direction op- 
posite to the cluster motion, is owed to the induced "wake" of 
background stars. If the star cluster is massive enough, the re- 
sulting deceleration can in principle be large enough to cause 
the cluster to spiral into the GC. 

Kim & Morris (2003; Paper I hereafter) performed numeri- 
cal simulations of dynamical friction suffered by star clusters 
near the Galactic center, and found that dynamical friction can 
indeed bring star clusters formed outside the central parsec 
into the central parsec. However, this is only possible if the 
cluster is either very massive (^ 10'' M0) or is formed near the 
central parsec (< 5 pc). In both cases, the cluster should have 
an initally very dense core (~ IO^MqPC"^). These extreme 
requirements make the dynamical friction scenario rather im- 
plausible. Clusters with smaller masses, larger galactocentric 
radii, and/or smaller core densities either 1) completely evap- 
orate before reaching the central parsec due to increasingly 
strong tidal forces that the cluster must endure during the in- 
ward migration, or 2) do not reach the central parsec within 
the lifetime of young, massive stars. 

McMillan & Portegies Zwart (2003) presented semiana- 
lytic calculations of the inspiral of star clusters near the GC 
to study the parameter space in which the cluster can reach 
the central parsec of the Galaxy within a few million years. 
While they performed simplified semianalytic calculations to 
explore a wide range of parameter space. Paper I implemented 
a numerical treatment that models both the cluster and the in- 
ner part of the Galaxy with a large number of particles to ac- 
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curately model the final dissolution phase of the cluster. The 
amount of mass deposited in the central parsec can be esti- 
mated more accurately in this way, and this estimation was 
the main goal of Paper I. 

Recently, Hansen & Milosavljevic (2003) suggested that 
if a cluster formed outside the central parsec harbors an 
intermediate-mass (lO^-IO'^M©) blackhole (IMBH),'* the dy- 
namical friction scenario may work with more plausible initial 
cluster conditions. The idea is that the deep potential well in- 
duced by the IMBH at the center of the cluster may be able 
to keep the cluster core intact against tidal disruption for a 
longer time during the inward migration by dynamical friction 
so that the cluster core can reach the central parsec. However, 
it is difficult to determine with the (semi)analytic approach of 
Hansen & Milosavljevic (2003) 1) how long the IMBH can 
hold the cluster core intact, and 2) how much stellar mass 
would ultimately be deposited into the central parsec after 
cluster dissolution. The present study carries out numerical 
simulations of dynamical friction on GC star clusters having 
an IMBH using the same numerical method as in Paper I in 
order to answer these questions. 

2. MODELS 

We use the same numerical models as in Paper I, except that 
the clusters in the present study harbor an IMBH at the center 
Here we briefly describe our models, and readers are referred 
to Paper I for a detailed model description. 

We use a parallelized N-body/SPH (Smoothed Particle Hy- 
drodynamics) code named GADGET (Springel, Yoshida & 
White 2001), which computes gravitational forces with a hi- 
erarchical tree algorithm and represents fluids by means of 
SPH. We implement only the gravitational part of the code. 

To model the potential of the central region of the Galaxy, 
we adopt a truncated, softened, spherical, power-law density 
profile: 
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This is a density model from Genzel et al. (1996) with an 
added exponential truncation, scaled by R,runc = 15 pc. The 
outer boundary of the Galaxy is set to be 25 pc. The net angu- 
lar momentum of the Galaxy is assumed to be zero. We use 
8.4 X lO'' particles to model the Galaxy, and the total Galaxy 
mass represented in our model is 4.3 x lO^M©.^ Thus the 
mass represented by a single Galaxy particle is ~ 5OM0. 

For the initial stellar density profile of the cluster, we adopt 
the Plummer density profile, 
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Clusters in our simulations initially have a total cluster mass 
Mci of 10^ Mq, galactocentric radius Rof 5 pc, core radius 
of 0.13 pc, and tidal radius r, of 0.77 pc (thus r,/rc = 6). The 
number of particles for the cluster is 10^, making the mass 
of each cluster particle ^ IOMq (this number is not exactly 
IOM0 because includes the mass of the IMBH, which is 
represented by a single particle). 

The Plummer density profile we adopted above is only for 
the stellar component, and we place an IMBH at the center of 

^ We do not address in the present study how such an IMBH would be 
formed. 

' Mgaiaxv for Models 2 & 3 in Table 1 of Paper I are erroneous and should 
read 4.3 X IO'^Mq and 6.0 X IO'Mq, respectively. 



TABLE 1 
Simulation Parameters 





^'inil 


MlMBH 


Pc 


rin 


Simulation 


(Vcrc) 


(Mo) 


(M0PC-3) 


(pc) 


1 


0.5 





lL4x 10' 


0.00 


2 


0.5 


2 X iC 


5.5 x 


0.06 


3 


0.5 


1 X 10" 


7.2 X lO*" 


0.05 


4 


0.5 


3 X 10' 


8.7 X 10'' 


0.04 


5 


0.2 


1 X Iff* 


7.2 X 106 


0.05 



this profile. However, an equilibrium configuration may not 
always be achievable when an IMBH is added to this profile 
(a part of the distribution function f{E) may become negative 
if the IMBH mass added is too large; see Binney & Tremaine 
1987 for details). It turns out that for the mass range of the 
IMBHs we consider in the present study, equilibrium is not 
available without modifying the adopted stellar density pro- 
files. We find that equilibrium is achievable when stars inside 
a certain small radius from the cluster center, r,„, are removed 
(once the simulation begins, the initial void of stars near the 
IMBH will be quickly filled by stars that are deflected into 
the void by two-body interactions, and the cluster will try to 
reach the new equilibrium). The values of r,„ that we choose 
are the minimum radii that allow equilibrium for the system, 
and are shown in Table ^ Each simulation has a different 
initial central density ipc) because of different r,-,, values. 

The initial orbital eccentricity is expressed in terms of the 
initial orbital velocity, v,„,v, relative to the circular velocity at 
a given R, Vcu-c, and our clusters initially have a tangential ve- 
locity only. Our simulations have v,„„ / VcVc of either 0.5 or 0.2. 
Such significantly non-circular initial orbits might result from 
cloud-cloud collisions, which is one of the possible mecha- 
nisms for forming a relatively massive cluster near the central 
parsec. Such highly eccentric initial orbits will make the dy- 
namical friction more effective, as the cluster will experience 
its largest tidal forces near its periapse position. Since this as- 
sumption favors the process we are investigating, it strength- 
ens our conclusion in what follows that the process is unlikely 
to be operating in the GC. 

As discussed in Paper I, clusters with > 10^ Mq would ini- 
tially contain more than a thousand YMSs, which is more than 
an order of magnitude larger than currently observed in the 
central parsec (there are very few YMSs outside this region 
except for two very young star clusters at least 30 pc from the 
GC with a mass of order of 10"*Mq, the Arches cluster and 
the Quintuplet cluster''). On the other hand, if Md is less than 
1O'*M0, the number of YMSs in the cluster would be only a 
fraction of what is currently observed in the central parsec, 
and one would need successive inspirals of several clusters to 
explain the observed number of YMSs. Thus we choose Md 
= IQr' Mq for our clusters. 

3. RESULTS 

Figure^compares the evolution of the radial distribution of 
cluster particles of Simulation 1 (Simulation 8 of Paper I) and 
Simulation 2, which have the same initial conditions except 
that the former does not have an IMBH. The orbital evolu- 
tion shows a periodic oscillation simply because the clusters 
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Fig. 1. — Grey scale map showing the temporal evolution of the radial 
distribution of cluster particles for Simulations 1 and 2. The IMBH is not 
present in Simulation 1 but is included in Simulation 2, though it is not shown 
in the plot. The grey scale bars next to the map represent the logarithmic 
scales of the density in units of Mqpc"'. The horizontal dashed line shows 
the location of R = 1 pc. In Simulation 2, the core (relatively bright, distinct 
spots) continues to sprial in below i? = 1 pc. 
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Fig. 2. — Evolution of galacocentric radius, R. of the IMBH (solid line) 
and the mass of the cluster stars that are located in the central parsec of the 
Galaxy, Mip^ {dashed line). 



initially have elliptical orbits. This plot dramatically shows 
the role of the IMBH in making the cluster core more stable 
against tidal disruption and prolonging its inward migration 
to smaller galactocentric radius. While the whole cluster of 
Simulation 1 completely disrupts when reaching R=l pc and 
is spread over the region between 0.8 and 3.5 pc, in the case of 
Simulation 2, where the cluster harbors an IMBH at its center, 
only the halo of the cluster disrupts outside the central parsec 
and the core continues to spiral in to the central parsec. Al- 
though the location of the IMBH is not plotted in this Figure, 
we find that the IMBH and the core are tightly bound to each 
other until the core is dissolved inside the central parsec at 
-3.5 Myr 

Paper I showed that without the IMBH, a lO^M© cluster 
must have pc of at least — lO'^Mopc"'' to deliver its core to 
the central parsec within the lifetime of the YMSs, but Simu- 
lation 2 lowers this requirement by almost two orders of mag- 
nitude. The initial pc of Simulation 2, 5.5 x lO''M0pc"-', can 
often be observed in the cores of dense globular clusters, and 



is close to the density estimated for the early phase of the 
Arches cluster (Kim et al. 2000). 

Figure |2j) shows the location of the IMBH and the mass 
of the cluster stars (not including the IMBH) located inside 
the central parsec, Mipc, as a function of time for Simula- 
tion 2. Mipc increases as the IMBH enters the central par- 
sec, and reaches its asymptotic value when the cluster core is 
completely dissolved. The stellar mass deposited in the cen- 
tral parsec is 8 x lO^'Mo for Simulation 2, which is ap- 
proximately three times larger than that of Simulation 1 (Fig- 
ure|2t). Clusters formed near the GC are expected to undergo 
very rapid dynamical mass segregation during the first few 
10^ yr (Kim, Morris, & Lee 1999). Thus by the time the core 
sinks to the central parsec and dissolves therein, the core will 
be largely dominated by the most massive stars. If we assume 
that the typical mass of a YMS is 50 M0 and the number of 
central parsec stars presumed to be YMSs is 50, then the to- 
tal mass of YMSs in the central parsec is 2.5 x 10^ Mq. This 
value is sufficiently less than Mipc at the end of the Simula- 
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tion that, allowing for a reasonable initial mass function, we 
expect that almost all YMSs formed in the cluster of Simula- 
tion 2 are transported into the central parsec. 

Mipc at the end of Simulation 1 is ~ 2.5 x IO-'Mq, which 
equals the estimated mass of YMSs in the central parsec. 
But only a few stars of Simulation 1 are located inside R = 
0.8 pc, while most of the observed YMSs are located inside 
R = 0.8 pc. The presence of the IMBH in Simulation 2 not 
only increases Mipc, but also transports the core stars deeper 
into the GC, as shown in Figure^ 

It may appear that Simulation 2 is a good scenario to ex- 
plain the presence of YMSs in the central parsec, but there 
is nonetheless a rather serious problem with it: the IMBH in 
Simulation 2 is too massive compared to the total cluster mass 
(20 % of the cluster mass). The typical core mass in a clus- 
ter is ^ 1 % of the cluster mass, so even if all the stars in 
the cluster core turn into an IMBH, it is still far smaller than 
the IMBH of Simulation 2. Furthermore, Portegies Zwart & 
McMillan (2002) estimate that an object resulting from the 
runaway growth through stellar collisions can grow only up 
to ^ 0.1 % of the cluster mass. The IMBH could grow fur- 
ther by accreting remnant gaseous material left over from the 
cluster formation before the gas gets blown away by massive 
stars, but it seems highly unlikely that there would be enough 
gaseous material for the IMBH mass to grow by a factor of 10 
to 100 by accretion in just a few million years. 

Clusters having IMBHs smaller than that of Simulation 2 
would be more realistic, but then the stabilization of the clus- 
ter core by the IMBH would be less effective. By comparing 
Simulations 2, 3, and 4 (Fig.lJJ), c, & d), which have the same 
initial conditions except for the IMBH mass, we find that the 
minimum IMBH mass fraction required to bring the core of 
the 10^ M0 cluster to the central parsec is 10 % (-- 10"* Mq). 
Mipc at the end of Simulation 4 is almost the same as that of 
Simulation 1, where there is no IMBH. It takes more than 
6 Myr for the core of Simulation 3 to migrate to the central 
parsec and dump the stars there, and this time is close to the 
upper limit of the estimated age of the YMSs in the central 
parsec. 

A more eccentric initial cluster orbit could be more effec- 
tive for the inward migration of the cluster core because the 
cluster would experience higher Galactic stellar densities at 
earlier times. However, Simulation 5 (Figure|2^), which has 
Vinit/vcirc = 5, shows that Mipc is not sensitively dependent on 
the initial eccentricity of the cluster orbit. 

We also performed a couple of simulations with an initial 
density profile from the King model (King 1966) with a con- 
centration parameter Wq = 9, which represents a highly con- 
centrated profile for a cluster, but found that the results are 
nearly identical to the simulations with the Plummer initial 
profiles. This shows that our results are probably not sensi- 



tive to the choice of the initial density profile. Initial R values 
other than 5 pc were not tried in the present study, but as found 
in Paper I, survival of the cluster depends on the initial pc, but 
not sensitively on the initial R. 

After trying and considering various different initial con- 
ditions for the cluster, we conclude that the presence of an 
IMBH that occupies less than ~ 10 % of the total cluster does 
not greatly increase the amount of stars deposited in the cen- 
tral parsec. 

The simulation method adopted in the present study prop- 
erly describes dynamical friction, but may not accurately fol- 
low internal dynamical evolution of star clusters due to the 
use of a rather large softening parameter (see Paper I), as well 
as to the uniformity of the model stellar masses, which cannot 
therefore reproduce mass segregation in the cluster. The two- 
body relaxation times at the core of our clusters are of order of 
10^ yr, so some cluster cores may experience significant dy- 
namical evolution during our simulation periods. In star clus- 
ters with a black hole of a considerable mass, the core density 
gradually decreases by the loss of stars to the black hole and 
the expansion of the core (Rauch 1999, Freitag & Benz 2002). 
This density decrease is probably underestimated in our sim- 
ulations, thus the actual initial core densities and black hole 
mass fraction required to explain the YMSs in the central par- 
sec are probably even higher than implied by our results. 

We note that from high spatial resolution millimeter line 
observations of the inner four parsecs of the GC, Christo- 
pher et al. (2004) found that some of the molecular clumps 
in the circumnuclear disk (1.2 > 7?/l pc > 2.5) appear to be 
dense enough to surpass the Roche limit in the central par- 
sec. The virial densities of their clumps range between 10^ 
and lO^Hcm"^, and the virial masses range between 10^ and 
10^ M0. Simulations in Paper I suggest that a cluster formed 
from the heaviest clump in the Christopher et al. observations 
(an inferred ^ 10^ M0) may be able to spiral into the cen- 
tral parsec within the lifetime of the YMSs, even without an 
IMBH. 
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